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1. Introduction 



In recent years a lot of work has been done aimed at understanding the physics 
of hadronic interactions in a dense (nuclear) environment. The questions ad- 
dressed in this context range form DIS at low x to high energy scattering 
on nuclear targets to heavy ion collisions. Experimental motivations come 
from the low x data measured at the DESY collider and the Fermilab Teva- 
tron and, of course, heavy ion collision experiments at RHIC and LHC An 
additional theoretical motivation is that the dense partonic environment pro- 
vides a link between straight perturbative QCD of factorization theorems and 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution, and the still 
poorly understood soft QCD physics. The main idea is that high partonic 
density might provide an effective infrared cutoff, so that the genuine infrared 
strongly interacting region becomes relatively unimportant. The contribution 
of the infrared region to various observables which are not infrared safe in the 
usual perturbative sense, then may be suppressed by powers of A//^, where 



1 



A is the infrared scale and /i is the scale roughly measuring the strength of 
the gluon fields in the dense partonic system, fi oc gA. If fi is large enough, 
the coupling constant governing important interactions is small, and many 
non infrared safe quantities may become in principle calculable. Although we 
are still far from the quantitative understanding of this interesting regime, 
progress is being continuosly made in different approachesErtj. 

In this paper our aim is to give a simple physical description and provide a 
transparent calculational framework to what one may call eikonal propagation 
approach. The eikonal picture views a high energy hadronic reaction as a 
projectile, consisting of a (usually small) number of partons impinging on a 
large target. The target could be either a nucleus in its rest frame, or a high 
energy hadron which has already developed many partons in its wave function. 
The interaction between the projectile partons and the target fields is assumed 
to be eikonal: that is the projectile partons propagate through the target 
without changing their transverse position but picking up an eikonal phase. 
As a result of the phase differences accumulated by the partons, the initial 
wave function decoheres and a different hadronic final state is produced. The 
target in this picture is described by an ensemble of gluon field configurations, 
and the averaging over this ensemble determines the cross section for various 
processes. 

Eikonal physics underlies recent approaches towards calculating the produced 
particle multiplicities in high energy hadron-nucleusi and nucleus-nucleusEa 
collisions. It is applied in the calculations of radiative jet energy loss in hot and 
cold nuclear matteJi as well as in many studies of low x DIS. A particularly 
lucid discussion with many explicit calculations is given by Buchmiiller et 
al.i. Moreover, as discussed in Ref.0, eikonal propagation is closely linked to 
the low X evolution of Refs.Bi. On the other hand, the eikonal approach has 
been critisized ini, since it does not take properly into account production 
of diffractive intermediate states, which should be equally important at high 
energyS. One of the main motivations for our work is indeed to set up a simple 
framework of eikonal propagation which can be extended to include diffractive 
production. While the present paper does not include diffractive production, 
we view it as a step towards this goaS 

This paper is organized as follows: In Section 2 we give an explicit quantum 
mechanical description of the eikonal evolution and show how to calculate var- 
ious cross sections from the knowledge of the outgoing wave function of the 
projectile that has undergone eikonal scattering. We explain how to perform 
averages over the gluon fields of the target. It turns out that the correct av- 
eraging procedure amounts to averaging on the amplitude level for the elastic 
and inelastic cross section, but on the cross section level for the diffractive 
one. 
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As a first illustration of our calculational scheme, we derive in Section 3 the 
gluon radiation in a hadron-nucleus scattering. This calculation reproduces 
the results of Ref.i for q — A and qq — A systems. 

In Section 4, we discuss Deeply Inelastic Scattering on a large target. Our cal- 
culations are consistent with the results of Buchmiiller et al.S, and provide the 
corrections to them. The most interesting quantity here is the diffractive 
cross section. We explain why physically the corrections in this case are 
precisely what distinguishes between the elastic cross section calculated iniS 
and the diffractive, or large rapidity gap cross section that we calculate here. 
We calculate the diffractive cross section due to the qq as well as to the qqg 
component of the virtual photon. 

In Section 5 we consider the low x evolution in DIS. For the total photoab- 
sorption cross section our results immediately yield the operator form of the 
low X evolution equation of Ref.i. We derive corrections to Kovchegov's 
version of this equation and show that these corrections express the fact that 
the fundamental dipoles do not propagate through the target independently 
beyond the leading order in 1/N. We also derive an analogous evolution equa- 
tion for the diffractive cross section generalizing the results of Ref.0. Here we 
again obtain the operator form of the evolution, which is a new result. We 
show that in the large N limit this equation reduces to that of Ref.ill and we 
derive all-order 1/N corrections to this equation within our model for the av- 
eraging over the target fields. We point out that both the leading and all-order 
evolution is not an independent equation. Once the evolution of the total 
cross section is determined, the diffractive cross section is obtained by direct 
integration of the known function of rapidity. 

Finally in Section 6 we summarize our results. 



2. Quantum mechanics of eikonal propagation 



Consider an energetic hadronic projectile impinging on a large nuclear target. 
The projectile is characterised by a wave function, in which the relevant de- 
grees of freedom are the transverse positions and colour states of the partons, 

= ^{{au Xi}) \{ai, Xi}) . (2.1) 



The colour index ctj can belong to the fundamental, antifundamental or adjoint 
representation of the colour SU{N) group, corresponding to quark, antiquark 
or gluon in the wavefunction. In what follows we will consider wave functions 
with a small number of partons. 
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A. The projectile biased observables 



At high energy the propagation time through the target is short, and thus 
partons propagate independently of each other. For the same reason the trans- 
verse positions of the partons do not change during the propagation. The only 
effect of the propagation is that the wave function of each parton acquires an 
eikonal phase due to the interaction with the glue field of the target. Thus the 
projectile emerges form the interaction region with the wave function 



with - the gauge field in the target and T" - the generator of SU (N) in 
a representation corresponding to a given parton. Thus the relative phases 
between components of the wave function change, and the state that emerges 
after the target is no longer an eigenstate of the strong interaction Hamiltonian 
(as the incoming state is assumed to be) but rather a superposition thereof. 

Given this explicit law for the evolution of the projectile wave function, one can 
calculate various observables related to the scattering process. For example, 
up to the factor of total flux, the inelastic cross section is proportional to the 
probability that the outgoing state is different from the incoming one. This 
probability is given by the norm 



(2.2) 



Here S is the ^'-matrix, and the W^s are Wilson lines 




(2.3) 



)Proj 
inel 



5^ 



(2.4) 



where 



[l-\^in){^in\]\'^out) 




(2.5) 



and s is the ^'-matrix element, or overlap 



s — 



(2.6) 



Thus 




(2.7) 
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The probability P^nd for a given configuration of tlie target fields of course 
depends on the impact parameter, that is the position of the projectile in the 
transverse plain. To get the cross section one simply integrates the probability 
over the impact parameter 

<c? = jd%PK^{h). (2.8) 



Note that the probability Pf^°\ as defined above does not include processes 
for which the projectile remains intact but the target is excited or breaks into 
fragments that stay close to the original rapidity. By its very definition, Pf^^i 
is the probability that the projectile scatters inelastically. We refer to it as 
a projectile biased observable and we denote it by a superscript "proj". The 



formulation which removes this bias will be given in section 2 C 



More generally, the incoming state may have internal degrees of freedom which 
we collectively denote by 7. This can be spin or (perturbatively) global colour. 
In that case the inelastic probability is the probability that the outcoming state 
is perpendicular to any one of the states \l/m(7)- Technically this means that 
5^ is defined by subtracting form '^out its projection onto a subspace spanned 
by 7 



m 



l-El^-(7))(*m(7)l 



\^out) 



(2.9) 



For orthonormal states ^'m(7) the operator 

Pi = l-El^-(7))(^m(7)l 



(2.10) 



is the projector on the subspace orthogonal to the possible initial states of the 
projectile hadron. 

One can also calculate other, less inclusive observables. For example the num- 
ber of gluons with momentum k produced in the scattering event is given 
by 

iVprod(k) = {5^a\k)a{k)\5^). (2.11) 



In principle, the calculation of any observable in an inelastic process has to 
allow for the evolution of the state after the interaction to the time when 
the measurement is performed, that is t = 00. Thus the relevant observable 
should be calculated not in the state |5\E' >, but rather in [7(0, cxd)|5\E' >, 
where ?7(0, 00) is the free evolution operator from the time zero (immediatelly 
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after the interaction) to infinity. Being a unitary operator, it cancels against 
its conjugate in the calculation of the cross section Eq. ( |2.4| ), but it affects 
some less inclusive observables, such as the number of produced gluons ( |2.11D . 
In the quantities we calculate in this paper its contribution is always a pertur- 
bative correction of higher order than the accuracy to which we are working, 
and we do not consider this type of correction. As will become apparent from 
our calculations in the later sections, this does not mean that we neglect con- 
tributions which in the Feynman diagram language correspond to "emission 
after interaction". Those in the present framework arise from the "overlap" 
contributions, proportional to the overlap s of Eqs. ( |2.5| , |2.(j| ). 

In contrast to the inelastic cross section, the total cross section includes the 
probability for an elastic process, for which the outcoming state is the same 
as the incoming one up to a phase shift. In this case, one includes in the 
calculation all components of the outgoing wave function which differ from 
the incoming one, 

|5^tot) = l^out) - l^in), (2.12) 

and 

CTtot = / d^b (^out(&) - ^in(&)|^out(&) - ^in(&)) 

= J (fh[2 - (^out(&))|^in(&)) - (^in(fe))|^out(&))] • (2.13) 

The elastic cross section measures the probability of a process in which the 
product of the collision is in the same state as the initial one but with a 
nontrivial phase shift. This is defined through 

|5*el) = [(^inl^out > -l]|*in), (2.14) 

and 

air = J d'b{6^,i\6^,i) = J d'b\l-{^,^,{bmUb))\\ (2.15) 



Another interesting quantity is the diffractive cross section. This includes all 
elastic and inelastic processes for which the outgoing state is different from 
the incoming one and is in a colour singlet, 

l^^diff) = ninglctl^out) - l^in) (2.16) 

and 
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^diff = J d% {Vsinglct"^ out{b) - '^in{b)\Vsinglct^out{b) - ^i„,(6)) . (2.17) 

The operator V is the projection operator on the colour singlet states. 



B. The target averages 



The only property of the gluon field in the target assumed in the preceding 
discussion, is that it does not vary during the time of the interaction. This is 
justified since the time of the propagation of an energetic projectile through the 
target is shorter than the natural time scale on which the target fields vary. The 
target is, however, characterised by a wave function which has a spread in field 
space. The relative phases of different components of this wave function are 
not relevant in the eikonal approximation, but the probability distribution for 
different field strengths is relevant indeed. In other words, while each scattering 
event happens on a frozen field profile, in order to calculate the cross section 
one has to average over the field ensemble which represents the wave function 
of the target. 

In this paper we adopt the same ensemble averaging procedure as in Refs.i'@'i. 
The target is pictured as an assembly of sources of the colour field extended in 
both the longitudinal and transverse directions. The sources are independent 
of each other and thus completely decorrelated. The longitudinal and trans- 
verse coordinates of the sources are randomly distributed within the target's 
longitudinal and transverse extent. The field emitted by each source is small, 
and to leading order in the field strength, the projectile can exchange only one 
or two gluons with each one of the sources. Pictorially, the average of, say a 
product of two Wilson loops can be represented as in Fig.l. 



= 1 



Fig. 1. The target average of the product of two Wilson lines, denoted by 
thick lines. The fields A'^ in VF's are distributed with Gaussian weights local in the 
longitudinal coordinate according to ( 2.18| ). 



Put into formulae, this means that the target average of the two-point func- 
tion of the vector potential A"*", which enters all our calculations through 
Eqs. dU,!;!) is 
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((A+'^(x, x-)A^\y, y-))) = 5'^' 5{x- - y') i?(x, y, a;-) . (2.18) 

All final results will depend on one particular linear combination of these 
averages 

t;(x,y) = j rfx-^[5(x,x,x-) + fi(y,y,x-)-25(x,y,x-)]. (2.19) 

For simplicity of notation we will sometimes consider the field ensemble to be 
translationally invariant in the transverse plain within some radius i?, and will 
take all fields to vanish outside this radius. In this approximation, v depends 
only on the difference of its coordinates. 

With this averaging procedure, the target averages can be explicitly calculated. 
This is done for a variety of observables in Appendix C. The simplest averages 
we will need are those corresponding to "fundamental dipole" and "adjoint 
dipole" scattering amplitudes 

l((IY[w^^t(x)iy^(y)])) =exp[-Cpt;(x-y)] , (2.20) 

^ {{i:i\w^\y)W\^)\)) = e^v[-CAv{^-y)] , (2.21) 



Ar2 _ 1 



with the Wilson loops in the fundamental and adjoint representations respec- 
tively. 

In most of the paper we will not specify the form of 5, but sometimes it is 
useful to keep in mind the following simple formiil 

i?(x,y) =x-y/i2(a;-), (2.22) 



where /x is a slowly varying function of x^ . The independence of /x on x inside 
the radius of the target corresponds to the leading logarithmic approximation 
in the following sense. We are working in the gauge where the only nonvan- 
ishing component of the vector potential in the target is . In this gauge the 
vector potential is related to the colour electric field by the following relation 
which is local in x~: 

A+(x) = r dxiEi. (2.23) 

The leading logarithmic approximation assumes that the scale of the variation 
of the color electric field is much larger than the transverse distances at which 
the field is probed (transverse momentum ordering). In this approximation, 
we have 
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A+(x) = (x-xo) -E 



(2.24) 



and thus choosing xq 

of Ei. [It is in fact not essential to choose xq 



0, ( pA8|) with (^ 



follow for isotropic distribution 
0, since the dependence 



on the base point xq will drop out from all physical quantities. All physical 
observables, as we shall see depend only on the difference A'^ix) — A'^{y) of 
the vector potentials.] In this leading logarithmic approximation, the dipole 
scattering amplitude of eq.( p.20| ) is 



exp 



-Qi 



(2.25) 



with the saturation momentum 
q2 = ^CfJ dx-iJ,\x~) . 



(2.26) 



Thus to reiterate, the procedure of calculating any observable in the eikonal 
approximation is the following. For a given initial state, we calculate its eikonal 
propagation according to (|2.2|). Next we "project out" the subspace which is 
spanned by the possible initial states, according to (|2.9| ). Then we calculate the 
observable as in ( ^.11| ). Then we perform the averaging over the initial state 
if this is needed - for example the averaging over the initial polarization in 
the case of unpolarized scattering. And finally we average over the target field 
ensemble in the two gluon exchange approximation with the only nontrivial 
target average specified in (p.l8|) (and (|2.22|)). 



C. With due respect to the target 



In the eikonal scheme described above, one does not follow the details of the 
dynamics on the target side, since the target is modelled by a field distribution 
rather than being assigned a wave function of its own. As a result, the inelastic 
cross section (|2.8|) calculated in this way is underestimated since it accounts 
only for processes in which the projectile scatters inelastically. Other inelastic 
processes in the usual sense, in which e.g. the projectile remains intact but 
the target breaks up, are not included in af^g/. 
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Fig. 2. Schematic picture of a difFractive contribution to the qq scattering cross 
section which does contribute to the inelastic scattering cross section. 



In the same way, the cross section o"g[°"' in ( [^.15| ) accounts for the probabihty 
that the projectile undergoes elastic scattering. Thus its definition includes 
the probability of inelastic processes where only the target is excited. Such 
"diffractive" processes are correctly included in the diffractive cross section 
CTdiff in (|2.17|) , but in the above, projectile biased discussion, they are shifted 
from their rightful place in amei into (Tei- As we shall see in our explicit calcu- 
lations, these contributions are suppressed by powers of in the situation 
where the projectile is much smaller than the target. It is however unclear 
to us whether the same is true when projectile and target are roughly of the 



same size. 



In order to have the proper separation between elastic and inelastic contri- 
butions, we consider now the role of the target wave function. We observe 
first that the target field averaging described above amounts to calculating 
averages of the appropriate Wilson loops in the target wave function. Thus 
for example 



((Tr [w^\^)W^{y)\)) ^ (vI/*jTr [w^\^)W^{y)\ 



(2.27) 



where \Ef*„ is the initial state target wave function which in the gluon field 
basis has the form 



\'^\n) = nm)- 



(2.28) 



We can now parallel the discussion of subsection |2 A| paying proper attention 
to the target wave function effects. The initial wave function of the projec- 
tile-|-target system is 



(2.29) 



{0!i,Xi} 
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The outgoing wave function is 

^^/^^t=S^^n= ^i{(^^^Xi})Il^ix^)a,(^^i^'[A]\{(3i,X^}) ®\A) . 

(2.30) 

Although the incoming wave function is a product of the wave functions of 
the target and the projectile, the outgoing wave function is not factorizable, 
since the eikonal factors W depend both, on the coordinates of the projec- 
tile particles and on the target fields. The 5- matrix element, or overlap that 
determines the inelastic scattering probability (|2.4| ) reads now 

S = i^^nl^out) = {[{W)t])p, (2.31) 



where the incoming and outgoing wave functions (|2.29| ) and ( p.30| ) contain 
both the projectile and the target parts, rather than just the projectile parts 
as in (|2.1|) and (|2.2| ) . Here we denoted the Wilson loop factors arising in the S- 
matrix collectively by W, {)t denotes the averaging over the target fields (wave 
function) as in ( p.27| ), and ()p is the expectation value with respect to the 
incoming projectile wave function^. The same goes for the elastic probability 
of eg. ( |2.14| ) and the total probability eq.( |2.12| ). Thus we find 



p 

inel 
Pel 
Plot 

inel , el , tot 



\l-s\\ 

Pinel + Pel = — S 
(f'bPinel,el,tot{b) . 



(2.32) 

(2.33) 
(2.34) 

(2.35) 



The S'-matrix element s, given by ( p.31|) , already contains the target average. 

We therefore see that in order to calculate the true elastic and inelastic cross 
sections, the "target field averaging" has to be performed on the anmlitude 
level rather than on the cross section level as is sometimes thoughtE3. The 
difference between these averaging procedures has a clear physical meaning: 
they differ in the treatment of the diffractive states close to the rapidity of the 
target (see Fig. ^. 

Not all quantities are affected by this more careful treatment of the target 
wave function in calculating the overlap ( |2.31| ). The total cross section is not 

^Here, in order to project out the projectile states with p^, we proceed in exactly 
the same way as explained in eqs. (|2j,|2l^). We thus assume that the only degrees 



of freedom that characterise the target are the vector potentials, and that there is 
therefore no need to perform an analogous projection on the target wave function. 
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sensitive to this, since it is linear in s and therefore contains in any case only 
one averaging over the target fields. Physically this is because it does not care 
about whether we assign the target side diffractive contribution to the elastic 
or the inelastic probability. The diffractive cross section is not affected since 
we defined diffractive processes as those in which the projectile wave function 
emerges from the interaction region in a colour singlet state. This definition 
does not refer to the final state of the target, and thus the target averaging 
should only be performed once, on the cross section level. The number of 
produced gluons is also unaffected if we are interested in gluons at rapidities 
far from those of the target. 

In the rest of this paper, we apply the eikonal approximation to a variety of 
simple cases. Whenever we are interested in the cross sections, we will calculate 
them according to both definitions: the one referring only to the projectile as 
described in subsection P A| and the one including the proper treatment of the 
target wave function, as described above. The projectile biased quantities we 
will denote by the superscript "proj". We calculate both quantities in order to 
illustrate the differences between them. We stress, however, that the physical 
inclusive cross sections are given by the expressions eq.( |2.35|) . In the eikonal 



approximation, the proper calculation of the elastic and inelastic cross section 
involves the target field averages on the amplitude level. 



3. Gluon production 



As a first illustration of our formalism, we present a toy calculation of the 
number of gluons produced in the high energy scattering of a quark off a 
nuclear target. 



A. Gluon production in q+A 



The incoming state in this case is a single quark idealized as the Fock state 
la) of the bare quark, supplemented by the coherent state of quasireal gluons 
which build up the Weizsacker- Williams field /(x). This can be written as a 
coherent state of a gluon field - the result of evolving the bare quark state by 
free evolution from time — cxd to the moment of the interaction 

U- \a) = exp J rfx|/(x)|2 + z J d^d^fi^) [a,(x, + 4(x, 0] |«) 

^\a) + |f/xf/e/(x)T^^|/3;6(x,0). (3.1) 
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Here, U- denotes the unitary free time evolution from — oo up to the tar- 
get, ad{'x,C,), a)^(x, ^) are the annihilation and creation operators for gluons 
of colour d at transverse positon x, and rapidity ^] and T'^ are generators of 
the fundamental representation of SU{N). Lorentz and spin indices are sup- 
pressed. We work in the light cone gauge appropriate to the projectile, that 
is A~ = 0. In this gauge the gluon field of the projectile is the Weizsacker- 
Williams field 

A*(x) oc^^(x-)/,(x), /^(a;^)oc^7^, (3.2) 



where x~ = is the light cone coordinate of the quark in the wave function. 
The integration over the rapidity of the gluon in the wave function (|3.1| ) goes 
over the gluon rapidities smaller than that of the quark. In the following we will 
suppress the rapidity label. We are interested in calculating the distribution 
of produced gluons per unit rapidity and since the wave function does not 
depend on rapidity, this distribution is obviously flat. 

The Weizsacker- Williams field is proportional to the QCD coupling and is 
considered to be small. The higher powers of |/(x)p are thus negligible and 
the approximation ( p.l|) applies. 

The interaction S of the projectile wavefunction ( ^.1| ) with the target results 
in the phase shifts described by eikonal Wilson lines, 

"^out = SU.\a) 

= W^<f,(0) |7) + / rfx/(x)T,%iyJ^(0X(x)|7;c(x)). (3.3) 



Here W^{0) and W'^{'x) are the Wilson lines in the fundamental and adjoint 
representations respectively, corresponding to the propagating quark at the 
transverse position Xy = and gluon at Xp = x. 

Gluon production being an inelastic process, we have to project out of the 
outgoing wave function the component lying in the subspace spanned by the 
incoming states with arbitrary colour orientation a. Thus the index 7 in the 
projection operator ( |2.10| ) has here the meaning of the quark colour index. 



(3.4) 



To lowest order in p only the no-gluon component of out contributes in the 
last (overlap) term. 



|5^„) = f/+ /rfx/(x 







<( 



-T'm^Ip^^) l7;c(x)) (3.5) 
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As illustrated by the diagrammatic representation of this equation, the struc- 
ture of the second term corresponds precisely to "emission after interaction" . 
By this we mean that the gluon is emitted from the initial quark state rotated 
by the eikonal factor W^, while there is no eikonal factor that accom- 
panies the gluon propagating through the target. Thus this term is the exact 
counterpart of the Feynman diagrams containing gluon emission from final 
state quark. Another point to note is that since the r.h.s. of ( |3.5|) has an ex- 
has to 



plicit factor /, to lowest order in / the evolution operator [/+ in 
be approximated by unity \cf eq. (|3.1|) ], as noted in the previous section. 



The number spectrum of produced gluons is obtained from ( |3.5| ) by calculating 
the expectation value of the number operator in the state 5\E'q,, averaged over 
the incoming colour index a. After some colour algebra (see appendix ^ for 
technical details) we obtain 



iVprod(k) = ^Y.i^^Mi^)(^d{k)\5^^) 



J rfxdye^''-(^-^)/(x)/(y) 

1 



((Tr [H^^t(x) 1^^(0)1 )) 



- 1 
1 



m - 1 



- 1 

((Tr [vr^t(y) l^^(O)])) 
((Tr [Py^t(y) W^^(x)l)) 



(3.6) 



Here, all the information about the scattering properties of the target are 
encoded in the target average of the product of adjoint Wilson loops which 
are determined by ( p.l8| ) and ( p.22|) . The explicit calculation of this as well 
as some more complicated Wilson loop averages is given in Appendix C. The 
result is 



1 



((Tr W^\y)W^{^) )K = exp [-C^ t;(x - y)] 



(3.7) 



or in the leading logarithmic approximation 



exp 



(x - y)^ Ca 
8 Cp 



Ql 



(3.8) 
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With the exphcit form of the Weizsacker-WiUiams field of the quark projectile 
in configuration space, 

Cf /(x) /(y) = {A,{x) A,(y))p,oj = ^ ^ , (3.9) 



the produced number spectrum (El6|) for g + A coincides with the expression 
given by Mueller and KovchegovO. Fourier transforming (|3.6| ), we see that it 
is nothing but the Gunion-Bertsch gluon radiation cross sectionE^I for a hard 
quark which receives a transverse momentum transfer k^, distributed with 
Gaussian probability distribution, 



iVprod(k) 



X I 1 + e 
SttCf 1 



2C 



k2 (k, - k)2 



(3.10) 



B. Gluon production in i^{qq) +A 



As a second example, we consider for the initial state wave function a quark- 
antiquark pair which may be viewed as a simplified representation of a hadron. 
We take the quark to be at the transverse position x and the antiquark at y. 
The state incoming into the interaction region is 



t/_ |(a,x)(a,y)) ^ |(a,x)(a,y)) + / c?z /(z - x) T„^(5s^ 



/(z-y)5„^r|^ |(/?,x)(Ay)(&,z)). (3.11) 



Here, the colour indices of the quark and antiquark are independent, i.e., the 
incoming quark and antiquark is not taken to be in a global colour singlet 
state. The reason for this is the following. Any realistic hadron wave function 
apart from a quark and an antiquark has a large number of soft gluons. Those 
are not the Weizsacker- Williams gluons present in eq. ( p.ll| ) but rather genuine 
soft nonperturbative glue. They do not scatter eikonally since their energy is 
small. One can reasonably expect that they are completely absorbed by the 
target. This absorption does not produce any gluons in the central rapidity 
region which we are attempting to calculate here. However the presence of 
these soft gluons serves to randomize the colour orientation of q and q. We 
thus model the hadronic state as having all colour orientations of individual 
quark and antiquark as equally probable. The implication of this is that the 
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projection operator which we employ to project out the component of the 
incoming hadron from the outgoing state, has independent summations over 
the quark and antiquark colour indices. 

The outgoing state is 

+ Jdz [/(z-x)T„%l^j;(x)W/^t(y)^A(^) 

-/(^-y)^|.<.(x)W^^^(y)W^i(z)] |(7,x)(7,y)(6,z)), (3.12) 



and after projecting out the initial state component 



5^, 



/ dz /(z - 



-/(z - y)7|. W^i;(x) W^^iy) W^,{z) 
-/(z- ' 



+f{^-y)T!;^w[,{^)w^{y) 



|(7,x)(7,y)(&,z)). (3, 



13) 



Here again we have kept only terms of the lowest order in The number 
of produced gluons is calculated in the state independently averaging over 
the initial state colour indices a and a. After some colour algebra we get 

iVprod(k) = ^ ^(5*„^|4(k) ad{k)\ 



Cf J dzdze^^^'-^'^ 

/(z-x)/(z-x) (l + e-^^(-^) 



^-Nv{z,-x) _ ^-Nv{z-x) 



+/(z - y) /(z - y) (l + e-^""^^-^^ - e-^''(^-y) - g-^^'^^-y) 



(3.14) 



In comparison with the number spectrum ( |3.(j| ), ( |3.14| ) is the independent 
sum of the radiation spectra for two quark-nucleus collisions, with quarks 
at transverse positions x and y respectively. [This can be seen by using 
t;(x - y) = ^^^Ql] In this sense, it appears sufficient to calculate the q-A 
gluon number spectrum ( p.6|) in order to determine the h-A number spectrum 
for "hadrons" h. However, there are at least two obvious limitations to this 
simplified treatment: first, if colour correlations between the projectile partons 
become important, the gluon emission off different quarks is not independent 
any more. This will be seen in our discussion of gluon production off a colour 



16 



singlet gg-Fock state in section Second, if the colour charge density in the 
projectile is large, then its WW field is not an incoherent superposition of the 
WW fields of its constituents. This in particular is important when both col- 
liding objects are nuclei, as in0. We plan to address this question in a future 
publication0. 



4. Small X DIS in the target rest frame 



In section |3 B| , we considered the propagation and accompanying gluon ra- 
diation of a quark- ant iquark pair in a nuclear target without requiring that 
this colour dipole is in a global colour singlet state. We argued that this is 
a simple ansatz for a meson-nucleus collision where the hard quark and anti- 
quark projectile components are not required to be in a singlet state since soft 
non-perturbative gluon fields are present as well. The situation is different for 
many problems in small x DIS. In the nuclear target rest frame, the process 
of interest here is the eikonal propagation of the hadronic Fock components 
of the virtual photon through the nuclear target field. Due to formation time 
arguments, these Fock states cannot develop a soft colour- carrying field prior 
to the interaction. The hadronic wavefunction of the virtual photon is hence 
perturbative, and the lowest lying qq and qqg Fock states of the incoming 
virtual photon are in a colour singlet. At not too low values of x = xq the 
photon wave function has only the qq component: 

= J c/2(x - y)dzij{x - y, z)-^Saa\a{x) , a(y), z) . (4.1) 

Here as previously x, y are transverse coordinates of the quark and antiquark 
respectively and z is the fraction of the longitudinal momentum carried by 
the quark. The explicit form of the wave function ip for both transverse and 
longitudinal photon is given in many places (e.g. Ref. and will not concern 
us here. We will also not indicate the z-dependence of the Fock states explicitly 
in the following, since the scattering amplitude of the qq pair in the eikonal 
approximation does not depend on z. 

The DIS cross section is determined by the total cross section of the qq pair 
on the target. The reason it is the total and not the inelastic probability that 
is important, is because even if only the phase of the qq state is changed in 
the scattering event, the coherence between the components of the 7* wave 
function is disturbed and thus what emerges from the interaction region is 
some hadronic state. The DIS cross section is thus given by 

a^'' = I dWyd^^(x - y, z)^*(x - y, z)PZix, y) . (4.2) 
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At lower values of x the virtual photon wave function develops a qqg Fock 
space component, due to the Weizsacker-Wiliams field of the quark and the 
antiquark. Thus the wave function of the qq pair in (|4.1| ) is substituted by the 
normalized state 

(5oa|a(x) , a(y)) = ^^(5os|tt(x) , a(y)) 



N VN 

In 1/x 



1= r^ddd'^[f{z-^)-fiz-y)] 

N Jlnl/xo J 

xr„%|(/?,x)(a,y)(a,z,0), (4.3) 



where in the leading logarithmic approximation we have taken the gluon emis- 
sion amplitude to be independent of the gluon rapidity ^. This wavefunction 
is normalized to unity via the factor 



:i + Cp/rf2zF2(x,y,z) ln(xo/a;)) 
Cf 



2 



J rf2zF2(x,y,z) ln(xo/a;). (4.4) 



Here and in what follows, we use the following notational shorthands for the 
Weizsacker- Williams fields, 

F ^ /(z - x) - /(z - y) . (4.5) 



The wave function ( [4.4|) is correct for Weizsacker- Williams photons with ra- 
pidity lower than that of both, the quark and the antiquark. The photons with 
intermediate rapidity can be emitted only by the fastest parton (be it q or q). 
In the leading logarithmic approximation, however, it is the soft photons that 
dominate the phase space and we only keep those in our analysis. 

Even though only P^^^ is relevant for inclusive DIS cross section, in the rest of 
this section we will discuss also the inelastic and diffractive cross sections of 
the qq scattering on the nuclear target. Those have distinct physical meaning 
and are useful for calculations of less inclusive quantities. Our starting point 
is the virtual photon component (|4.3| ). The interaction with the target evolves 
its wave function into 



^out = S ^i. 



N 
1 



N 



W^\y)W^{-4\_ |a(x),a(y)) 



W^\y)T-W^{^)\_ iy,^,(z)|a(x),a(y),6(z)), (4.6) 
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where we have suppressed the integration over the gluon rapidity. 



A. qq- contribution to photoah sorption cross sections 



We start with the discussion of the contributions from the lowest lying qq Fock 
state. Details of the medium-average as well as more differential cross sections 
are listed in the appendices. 



The total scattering probability, which determines the inclusive DIS cross 
section according to ([4.2|), reads 



-'tot 



((2 



1 

N 



Tr 



1 

N 



Tr 



-Cf t'(x-y) 



(4.7) 



For the inelastic photoabsorption cross section, we get 
1 



inel 



1 - — ((Tr [W^\^) W^^(y)J)) ((Tr [W^^{y) W^{^)\)) 

1 _ e-'^Cp v{x-y) 



(4i 



and 



Mnel 



{{l-±Tr[w^\^)W^iy)] Tr 

/V2 _ 1 1 



Ar2 



-Nv{x-y) 



(4.9) 



As explained in section |], the difference between these two expressions is due 
to the target side diffractive processes which are not included in ( |4.9| ) but are 
accounted for in ([4.8| ). Comparing ( [4.8| ) and ([4.9| ) we see that this difference 
is of the order 0(1/A^^). The diffractive probability reads 



pqg 

diff 



(((lTr[H^^(x)l^^t(y) 



1 



A^ 



((^ (l + W^a1(x)W^.1(y)) - ^Tr 

\W{y) iy^t(x 



Tr 

N 



N' 
1)) 



1)) 
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-Nv(yL-y) 



(4.10) 



The leading O {^j^^ contribution in ( |4.10|) is the result given by Buchmiiller 



et al. for the diffractive photoabsorption probabilityi. It is easy to see from 



eqs. ( [4. 8| , [4 .71 , ^.351) that this contribution is precisely equal to the elastic proba- 



bility Peu since Pei = Ptot — Pinei- This leading O (^j^^ result was later rederived 
in RefS and it lead Kovchegov and McLerran to conjecture that the diffrac- 
tive cross section can be calculated by the procedure of averaging over the 
target fields on the amplitude level. Physically, however, as noted in RefJ§ 
and explained above, such a procedure corresponds to (quasi)elastic scatter- 
ing where each target colour source (nucleon) in the target nucleus remains 
in a colour singlet state. In contrast, the diffractive production requires the 
outgoing projectile wavefunction to be in the colour singlet state but it can 
leave some target sources in the colour octet and thus the whole nucleus in 
an excited state (the overall colour neutrality of the target is of course still 
preserved). Our derivation makes it clear that for the diffractive cross section 
the averaging must be done on the level of the cross section. We also see that 
the higher order in 1/iV corrections have a clear physical meaning in this case 
- they distinguish between the diffractive and the (quasi)elastic production. 

As our previous discussion suggests the difference between the diffractive and 
elastic cross sections for the scattering of qq pair is due to precisely the same 
physical final states as the difference between Pjnei and -Pjne'i''- Indeed subtract- 
ing from the total cross section (|4.7| ) the projectile biased inelastic one ( |4.9|) , 
one obtains the diffractive cross section ( [4.10| ). For more complicated projec- 
tile configurations this is not necessarily the case, since the global color is not 
the only degree of freedom that characterizes the outgoing state. This is also 
not true beyond the eikonal approximation where the transverse positions of 
the partons in the projectile may change during the interaction time@. 



B. qqg- Contributions to photoabsorption cross sections 



The non-abelian Weizsacker- Williams field in the virtual photon wavefunction 
( |4.3| ) is perturbatively small. In the calculations of the last subsection, this was 
used to neglect all powers of F. Now, we discuss the 0(F^)-corrections to these 
photoabsorption cross sections. There are several motivations to do so. Though 
nominally suppressed by powers of the strong coupling constant, such radiative 
corrections are known to give important contributions to the leading twist 
photoabsorption cross sections in some regions of phase space. Moreover, the 
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gluon radiation off the gg-Fock state can be viewed as a leading contribution 
to the small x evolution equations. This we shall discuss in section ^. The 
0(F^)-corrections arise from the gg^f-components of the ingoing and outgoing 
7* wavefunctions (^), (0) and (^). 



The total scattering probability reads 



tot 



N 



+ ln(xo/x) {{\ ^Tr [w^\^) W^{y)] + ^Tr [w'^^y) W^{^) 



2N 
1 



-Tr 



■CFv(yL-y) 



= (2 - 2e" 
+2Cf \n{xo/x)F'^ (e" 



2N 



:Tr 



ly^^(z) W^^(x)l Tr \W^\y) W^{z 



W^^i^) W^{z)] Tr \w^^{z) W^{y)] \)) 



-CFi>(x-y) _ -f [i,(z-y)+t,(z-x)] + 5i^?;(x-y) 



(4.11) 



Here and in the following we use the shorthand 



F'^ |rf2zF^(x,y,z)[... 



(4.12) 



which stands for the operator acting by multiplication with the square of the 
WW field, and integration over its coordinate. Thus in this and the following 
expressions the coordinate of the gluon z is always integrated over. 

The inelastic photoabsorption probability reads 



p: 



1 



inel 



Ar2 



1-— ((Tr W^^{^)W^{y) ))((Tr W^^{y)W^{^) )) 



1 



+ Hxo/x)F' ^ ( (Tr [W^^ (y) W^{^)\ ) ) 



x{{\NTT\w^\^)W'"{y) 



Tr 



PF^t(x)l^^(z)l Tr [iy^t(2)|y^(y)l U) 



1 



+ H^o/x)F'^^{{Tt [W^\^)W^{y)\)) 



x((|iVTr [W^^y) 



Tr 



iy^t(z)iy^(x)l Tr [w^^^(y)W^^(z)l )) 



1 _ ^-2Cpv{^-y) 

+2C^ln(a;o/a;)F2e-^^"("-y) |^g-c^.(x-y) _ ^-f [.(z-y)+.(z-x)]+^.(x-y)^ _ ^4_^3^ 



Finally, the diffractive probability reads 
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p{gg+ggg) _ / / f ^ rp 
+Cf ln{xo/x)F'- 



1 



+ 



m - 1 



Tr 



1 



W^\y)W^{z)\ Tr \w^\z)W^{^ 



N{N^ - 1) 
1 



Tr 



X 



iV2 - 1 
1 



Tr 



W^\^)W^{z) Tr W^\z)W^{y) 



- 1) 



Tr 



W^{y)W^\^)]-l) )). (4.14) 



The explicit averaging of this expression over the target field is calculated in 
Appendix 0. The final result is 



^diff - l^Hv^ + 



Ar2 ' 



-ArD{x-y) 



-Cf ln(a;o/x)F2 
iV2 _2 



A^2(Ar2 _ 1) 



Ar2(iV2 - 1) ^ iV2(Ar2 _ 1)2 

iV2 -4 



2(Ar2 - 1) 



-Af[u(z-x)+i.(z-y)] 



Ar2(Ar2 _ 1) 



g-f^(x-y) 



/AT + 3 



4(iV2 _ 1) VA^+ 1 



-(Af+l)[i)(z-x)+i'(z-y)]+t;(x-y) 



' N -I 



-(Ar-l)[«(z-x)+i)(z-y)]-i)(x-y) 



(4.15) 



To leading order in 1/A^ this simplifies to 



P. 



diff 



1 _ g--*'{x-y) 
+CfF^ ln{xo/x) |- (l - e-^''(""y) 



1 _ g-f b(z-x)+^,(z-y)]^2| 



(4.16) 



For completeness, we also give the projectile biased inelastic probability 
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2 ^ ^ 



+ ^iy/,(z) W^,^,(y) iy,^/x) + /ace/fcd/) 



+2C^ln(xo/x)F2|i, 



L ^ -NvU-y) '^^Nv{z-yz) 

_}_ Nv{z-y) _ f-i _ AS] p-f [f(z-x)+i,(z-y)+i,(x-y)] 1 



(4.17) 



To leading order O (;^), the qq contribution to the projectile biased and true 
photoabsorption cross sections coincide. However, for the qqg configuration, 
the elastic and the diffractive cross sections do not agree even at the leading 
order in 1/A^. The elastic cross section can be immediately inferred from ( [4.111 ). 
To leading order in 1/A^ 

p{qq+qqg) = [l - g-f ^(^-y))^ 

+0{l/N^). (4.18) 



This has to be compared with ( |4.16|) . The difference in the O(F^) terms is 
due to the fact that as opposed to the qq case, the quantum state now is 
characterized by a new degree of freedom. This degree of freedom is the relative 
phase between the qq and the qqg components of the wave function. Thus the 
outgoing state can be a colour singlet which is not identical to the incoming 
one. 



5. Eikonal low x evolution 



The results of the previous section allow us to derive evolution equations for 
inelastic and diffractive DIS cross sections in the eikonal approximation. In 
the operator form, the equation for the evolution of the inelastic cross section 
is the same as that of BalitskyS. Our derivation is very simple and allows us 
to trace explicitly some of the effects that limit the range of validity of this 
equation. We also present the all order 1/A^-corrections to Kovchegov's version 
of this evolution equation. Those are obtained with the use of Fierz identities 
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as described in the appendices. The equation for the diffractive cross section 
as given here in the operator form is new. In the large N hmit it reduces to 
the equation discussed in Ref.0 and we present corrections to this form. 



A. The inelastic cross section 



Consider DIS at some low value of Xq. We work in the frame where the pho- 
ton wave function at Xq has only a qq component, but the target ensemble 
is characterised by large values of the field. There are two distinct physical 
situations that correspond to this. One is the scattering on a nuclear target 
at not very high energy. In this case the large target fields are due to large 
number of nucleons in the target rest frame, and no significant x evolution 
has taken place (xq is not very low). Another situation is that the scattering 
is on a single hadron and Xq is very small. In this case the large target fields 
have been generated due to the low x evolution which in our frame has been 
put into the target. There are some clear physical differences between the two 
situations, which we will discuss later, but for now we treat them both in the 
same framework. 



Thus at Xq the DIS cross section is given by ( [4.2|) . The first line of (|4.7| ) then 
expresses it in terms of the glue fields of the target, while the second line gives 
its explicit expression in terms of our ansatz for the target field averages. Now 
consider a lower value of x. It is convenient to use the frame where all the 
evolution between xq and x is occuring in the photon wave function. This 
means that to achieve higher interaction energy the photon is boosted, while 
the target remains the same. Thus the target field ensemble is unchanged, 
but the photon wave function at x has an extra qqg component according to 
( |4.3| ) and ( ^4.4| ). The inelastic scattering probability is now given according to 
(gj) but with P^m+qqg) of substituted for P^^". If we were to put the 

evolution in the target instead, the inelastic cross section would still be given 
by (|4.2| ) but with the average performed over the new target field ensemble. 
We can thus combine ( [4.7| ) and ( 4.11 ) into a differential equation governing 
the change of the target field averages with x (we define ^ = ln(a;o/a;)) 



[W'''t(x) W''{j)\)) = -F=({yTr [w"'»(x)W'^(y) 



-It. 

2 



W^\^)W^{7)\ Tr [w^^^(z)W^^(y)l)) . (5.1 



This equation is identical to that derived by BalitskyB. Taking now the av- 
erages over the target fields and denoting the "dipole scattering probability" 
by 



24 



D(x,y) = 1 -e-^^^^^'-y), 



(5.2) 



we have 



^D(x,y)=C^F2|l-D(x,y) 

-[1 - D(z,y)]l^[l - D(x,z)]^[l - D(x,y)]l^ . (5.3) 



To leading order in and assuming that D does not depend on the impact 
parameter (6 = (x + ?/)/2) so long as h is smaller than the target radius R, 
this reduces to the equation derived by KovchegovB. 

The extra terms in ( |5.3| ) furnish corrections to Kovchegov's form of the 
evolution equation. Interestingly enough the corrections arise as noninteger 
powers of [1 — D], and thus at small D correspond to infinite series. In the 
saturation regime, where D is close to unity the effect of these corrections may 
very well be large, and it would be quite interesting to explore these effects 
numerically, similarly to Ref.il. 

These corrections have a very simple meaning. As discussed by KovchegovS, 
the leading approximation in the present context is equivalent to assum- 
ing that the dipoles interact with the target independently. From our deriva- 
tion of the target averages in Appendix |C|, it is obvious that the corrections 
in (|5.3| ) arise precisely due to the breakdown of the independent propagation 
beyond the leading order in 

What is meant by "independent propagation"? Physically, as discussed in 
section ^, in the eikonal approximation the partons of the projectile propa- 
gate through the target field completely independently of each other. This 
statement follows from the eikonal kinematics and does not rely on the 
expansion. However the fields in the target are correlated with each other at 
different transverse coordinates. As a result the eikonal phases picked up by 
the partons are not independent when averaged over the target field ensem- 
ble. If the target fields are small locally at every longitudinal coordinate, so 
that our model for averaging is applicable, then in the leading order in 
all eikonal phase averages factorize into amplitudes corresponding to funda- 
mental dipoles. In this sense, the dipoles "propagate independently" only to 
leading order in 1/N. Moreover, it is crucial for this factorization that any 
given source in the target can exchange at most two gluons with the projec- 
tile. If the fields are locally not small it is possible that the propagation is still 
eikonal, but the factorization breaks down. 

We finally discuss the limitations of this approach. There are at least two 
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reasons to expect that the evolution equation derived here breaks down at 
very low values of x: 

Firstly, the derivation assumes that all the evolution up until xq occured in the 
target wave function. Further evolution is nothing but boosting of the target. 
In the boosted frame the gluon fields contract in the longitudinal direction and 
their amplitude grows, as discussed in Refs.0i. This does not mean that the 
effective width of the target shrinks, since the low momentum components of 
the glue field also become important. The effective width of the target then can 
remain approximately constant, but the fields inside it grow. Therefore if the 
target wave function has undergone a lengthy low x evolution, the target fields 
have significantly increased locally and the averaging procedure employed here 
and in Ref.B breaks down. Thus we expect (|5.3| ) to be a sensible description 
for targets which are large but "soft", like a large nucleus in its rest frame, 
but not for targets which are "hard" , like a very energetic proton. 

Secondly, if the evolution is followed to very small values of x, the projec- 
tile wave function itself becomes dense. Then any further evolution has to 
take into account that the Weizsacker- Williams field of the projectile becomes 
large. Further evolution then becomes nonlinear already on the projectile side, 
since we would not be able to neglect higher orders of the field F in all our 
calculations. The onset of this effect may be expected to occur at the scale 

ln(xo/x) ~ ^ , (5.4) 



when the qq and qqg contributions to the photoabsorption cross section become 
comparable. In the language of the second paper in Ref.Q, such higher 0(F2)- 
corrections are due to "Pomeron loops". The evolution equation then breaks 
down already in the operatorial form. Another way of saying it, is that the 
averages of the products of large number of Wilson lines is not described by 
a straightforward generalization of (|5.1|) given in Ref.0. It is an interesting 



question whether the evolution equation ofH accounts better for this physics. 
As discussed in0, the JKLW evolution is different form the BK evolution, 
when the fluctuation fields are large. This is precisely the case of a large 
projectile. However whether this difference is due to different physics content 
or to technical__difficulties with JKLW equation is not clear at the moment 
as discussed ii 



^ It has been recently suggestedcj that the apparent differences between the BK 
and JKLW equations are due to the incorrect derivation of the latter. However, 



since the discussion ofS does not point tQ_a mistake in the calculations oi^"^^ 
consider the question open as discussed inll3. 



we 
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B. Diffractive cross section 



We can repeat the same excercise for the diffractive cross section. Using 
eqs. ( |4.iq ) and ( |4.14| ) we get 

PM^,y) = -CfF' (((^Tr \W^{^) W^\y)] - 1 



X 



lTr[iy^(y)iy^t(x) 



1 



- 1 



Tr 



1 



X 



N{m - 1) 
1 



Tr 



- 1 



Tr 



1 



W^\^)W^{z) Tr W^\z)W^{y 



N{N^ - 1) 



Tr 



W^{y)W^\^)]-l])). (5.5) 



This is the operatorial form of the low x evolution equation for the diffractive 
photoabsorption cross section. 

A comment is in order about the precise meaning of Pdifi- We have defined Pdifr 
at Xo as the probability that the qq pair in the photon wave function remains in 
a colour singlet state after the interaction with the target. At lower value of x 
the wave function contains a quark- ant iquark-gluon component. The rapidity 
of this Weizsaker- Williams gluon in the leading logarithmic approximation is 
always smaller than the rapidity of the quark or antiquark but greater than 
the initial rapidity ln(l/xo). Again we require that the outgoing state is in 
the colour singlet. This requirement does not constrain the emission into the 
rapidity interval between the gluon and the quark- antiquark pair in the final 
state, but it does not allow emission of particles with rapidity smaller than 
that of the softest gluon, that is ln(l/a;o). Thus the final states allowed by our 
definition have rapidity gap on the target side of at least ln(l/a;o)- 

Our definition of Pdifr coincides therefore with that of Ref.0 and (|5.5|) is 
the operator form of the evolution equation for this quantity. It is also 
easy to see that in the leading order in this equation indeed coin- 

cides with the equation derived in Ref.l2il. The straightforward way is again 
to realize that in the leading order in the target field averages fac- 
torize. Thus the right hand side of (|5.5|) can be factorized into products 



Ij (h.c.))) (denoted inll by A^^) and 
)) (denoted in§ by A'o). Neglecting the 
explicit 1/A^^ corrections in the coefficients of (|5.5|) the result is identical to 



of PdifT = (((^Tr 
Ptot = (((l-^Tr 



Vr^(x) W^^y) 
W^^(x) W^^iy) 
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the differential form of eq.(ll) o: 



d 



Pdis{x, y) = CpF^ {-Pdis{x, y) - Pd\s{x, z) - Pdm{y, z) 

+Pdifr(a;, 2)Pdifr(^, y) - 2Pdifr(a;, z)Ptot{z, y) 
-2Ptot{x,z)PAm{z,y) + 2Ptot(a;, 2;)Ptot(^, 1/)} • 



di 



(5.6) 



We note however that within the region of vahdity of the present approxima- 
tion, (|5^ ) is not really a differential equation, but it rather determines Pdifr 
directly in terms of the dipole cross scattering probability D of ( |5.2|) . The 
point is that the right hand side of (|5.5|) in the leading order in can be 
written entirely in terms of D: 



^Pdifr(x, y) = CpF' {-D\x, y) + [D{x, z)D{z, y) - D{x, z) - D{z, y)f} . 

(5.7) 



Thus once -D(0 is found by solving (|5.3| ), the diffractive probability can be 
found by directly integrating the right hand side of (|5.7|). The reason ( |5?7| ) is 
consistent with ( ^.6|) is the following. The initial condition satisfied by Pdig as 
discussed inEllis 



At lower values of x the solution of (|5.71) of course will not satisfy this relation 
any longer. However the difference between Pdiff and at any value of x 
is small, of order ln(xo/x)P^. Thus in the right hand side of the differential 
equation ( |5.6| ) the two can be interchanged as long as ln(xo/a;) <^ 
As discussed in the previous subsection, this condition has to be satisfied 
anyway for the whole scheme to be consistent. For even lower values of x the 
corrections nonlinear in P^ (or "Pomeron loop" diagrams) become important, 
and those are not accounted for in the present approximation (as well as in 
the derivations of Refs. 

Thus we stress again, that in order to find the diffractive probability one does 
not have to solve an independent evolution equation, but rather integrate 
(|5.7| ) with the known right hand side. This statement remains true beyond the 
leading order in The complete set of corrections due to correlated 
dipole propagation is read off eq. ( |4.15|) , 



Pdifr(x = xq) = D{xof . 



(5.8) 



d_ 

d^ 



PMx,y) = CFF' 



Cf f (x-y) 



(iV2 _ 1)3 _ 1 

Ar2(Ar2 _ 1)2 



Aft)(x-y) 



A^2(Ar2 _ 1) 



+ 2e 
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-2 



JVt)(z-x) 



-Nv{z-y) 




Ar2 -4 



Af[t)(z-x)+i;(z-y)] 



Ar2(Ar2 - 1) 



+ e 



_2g-f b(z-x)+Kz-y)] 




Ar2(Ar2 _ 1) 



2 



f^{x-y) 



) 



4(A^2 _ 1) Va^+ 1 




(Af+l)[i){z-x)+-y(z-y)]+-u(x-y) 




(Ar-l)[»;(z-x)+t){z-y)]-i){x-y) 



Here, y(x 



y) is related to the total dipole cross section D(x — y) via (|5.2|). 



6. Summary 

In this paper we have discussed in a very explicit and simple way the physics 
of high energy processes in the eikonal approximation. This provides a simple 
unifying framework for a variety of calculations existing in the literature. In 
particular we considered processes where a small projectile (a hadron or a 
virtual photon) scatters on a large target (a nucleus). Our approach then al- 
lowed us to derive the inclusive one particle spectrum as well as total, inelastic 
and diffractive cross sections. In particular the low x evolution equation for 
the DIS cross sectionB follows immediately from this calculation. The various 
sources of the breakdown of this evolution equation at very low x are also seen 
clearly in our derivation. We have also derived the operatorial form for the 
evolution of the diffractive cross section and have shown that in the large 
limit it reduces to the equation of Ref.0. We have clarified that this evolution 
is not independent of the total photoabsorption cross section, and indeed the 
diffractive cross section is determined by a direct integration once the total 
dipole scattering probability (at fixed impact parameter) is known. 

The averaging over the target nucleus wave function has been performed in 
the approximation where the target is considered to be a dilute system of 
weak sources of the gluon field. We did however calculate explicitly all-order 
corrections in this model. This allowed us to provide all-order corrections 
to the various cross sections as well as to both evolution equations. These 
corrections stem from the fact that fundamental dipoles do not propagate 
through the target independently from each other. These corrections have not 
appeared in the literature so far. 

We hope that this work is useful in providing a simple explicit unifying frame- 
work for the eikonal physics. The method discussed here allows simple calcu- 
tions of many other observables at high energy. What we consider to be most 
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interesting are less inclusive observables giving detailed information about the 
structure of the final states. It is for example straightforward to determine evo- 
lution of such quantities as the inclusive one particle distributions and cross 
sections for diffractive vector meson pro duct ionli3. 

Acknowledgement: We thank Yuri Kovchegov and Sangyong Jeon for in- 
teresting discussions and for pointing out typos in an earlier version of this 
manuscript. 



A. Photoabsorption cross sections 



In this appendix, we give full expressions for the differential inelastic (c = inel), 
diffractive (c = diff) and total (c = tot) photoabsorption cross sections. We 
consider contributions from the qq- and qqg- Fock components of the virtual 
photon. We also give expressions for the differential cross sections with respect 
to the transverse momenta pi, P2 and k of the quark, antiquark and gluon 
respectively. The calculation of the medium average {{■■■)) will be discussed 
in Appendix |C[ 

The wavefunctions S'^c are constructed as discussed in section ^ 

The ^-matrix element, or overlap (p^.6D at fixed transverse coordinates reads 



(VI/ 



(gg+ggg) , 



x,y)|^ 



out 



x,y)) = (r-^Tr[l^^t(y) w^{^ 



+ lln(xo/x)|ci2z|/(z-x)-/(z-y)p 

[T'^W^i^)T'W^\y)] W^,iz))). (Al) 



xTr 



Treating higher powers of the WW fields / in the normalization factor Af as 



perturbatively small contributions which can be neglected, and using (|B.5|) 
one finds 



x,y)|^ 



1 



out 



(A.2) 



2N 



^-Hxo/x) J rf2z|/(z-x)-/(z-y) 



NTt 



For the diffractive cross section ( p.l7|) , we use the projection operator on the 
colour singlet states 
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(A.3) 

(A.4) 



We start with the contribution from the qq component. The contribution to 
the total cross section reads 



(A.5) 



For the inelastic cross section we have 



(5*?lei(x,y)|<J*fl(x,y)) 



W^\^) W^{y)])){{^Tr [w^\y) W^^(x)])) , (A.6) 



while the projectile biased inelastic probability is 

(i-I'?»r'(x,y)l**?«er(x,y)) = 



ly^t(x) W^{y)] Tr [w'^^y) W^{x)])) . 



Finally the diffractive probability is 

(5^fff(x,y)|5^^Sr(x,y)) 
= K(x) w^^^y)] - 1 



(A.7) 



— Tr 



W^{y)W^\^)]-l^)). (A.8) 



For the same expressions including the qqg component, we find: 

{S^lT'''\^,y)\6^[T'''\^,y)) = 



2N 



((2 Tr 



X jTr [W^^{z) ly^(x)] Tr [w^^^(y) W{z)] - NTr [w^^^(y) ly^(x)] 

+Tr [W^\x) W^(z)] Tr [w^\z) W^(y)] - NTr [w^\yi) W^(y)] \)) , (A.9) 



{Mlir''\^,y)miir"{^,y)) 



(qq+qqg) , 
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+ Hxo/x)F' ^ ( (Tr [W^^ (y) W^{^)] ) ) 



x((|iVTr [iy^t(x)iy^(y) 



Tr 



1 



Tr 



W^\z)W^{^)]TT\w^\y)W^{z)] )), (A.IO) 



(gq+qgg) , 



x,y)) 



((f^Tr [ly^(x) ly^^(y) 



-Tr 



+Cf ln(xo/x)F2 <! - ( ^Tr [l^^(x) W^\y 



+ 



1 



A^2 _ 1 



Tr 



W^\y)W' 



Tr 



N{m - 1) 
1 



Tr 



X 



- 1 

1 



Tr 



Tr 



N{m - 1) 



Tr 



(A.ll) 



The virtual photon wave function in the approximation ( [4 .51) contains at most 
one quark, one antiquark and one gluon. Thus the various differential cross 
sections are obtained form the appropriate number of particles by integration 
over the impact parameter. For example 



dpi dp2 



rf^x (x, y) |aT (p0a,(pi)4(p2)s-p2) |5vl>f (x, y)) 
J dzd^^ d V d^y e^P^ (^-*)+^p2 (y-y) 

xV^(x - y, z)^*(x - y, 2;)Pe(x, y, x, y) , (A.12) 



and 



da™ 



(ipi (ip2 (ik 



j d^ dy dz dx dy dz e^P^ (''-*)+*p2 (y-y)+»k (z-z) 
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x(5v^r(x,y,z)|at(p0a,(pi)at(p2) 

xa,~p2)at(k)a(k)|5^f^(x,y,z)) 
dzd^^ d^z d^x d^y d^z e'^^ (y~y)+^^ ^--) 

xV^(x - y, z, z)iIj*{-k - y, z, z)Pc{5t, y, z, x, y, z) , 



(A.13) 



where Og, aq and a are annihilation operators of the quark, antiquark and 
gluon respectively. The wave function ipiii. — y,z) is defined in ( [4 .11 ) and 
V'(x — y, z, z) can be read off ( [4 .11) and (|4.3|) . Pc are amplitudes given in the 
following in terms of the target gluon field averages. We give the expressions 
for the quark-antiquark-gluon probability Pc(x, y, z, x, y, z), from which one 
can obtain the quark-gluon probability as 



Pc(x, y, X, y) = y d^zPc(x, y, z, x, y, z) 



(A. 14) 



We find 



Ptot(x,y,z,x,y,z) 



{{MM ( -Tr 



W^^t(x) w^{y) W^^y) Vr^(x) 



iTr 



+ 1 



+ —FF 

N 



l^-(z)l^/,(z)Tr 



-Wi(z) Tr 



-W^a1(z)Tr 



+ 



AT^ 



(A.15) 



Pinci(x,y,z,x,y,z) = (( 



N 



Tr 



W^\^) W^iy) W^\y) l^^(x) 



1 ^ 



N 



Tr 



W^\y) ly^(x) 



i?2r|n 

Ar2 



W^\y.)W^{y) \ Tr [T"iy^(x)T^iy^^(y)l iyi(z 



iV2- 



P^Tr [iy^(y)T''iy^t(x)T"l Wi{z)llx\w^\y)W^ {^^)) , (A.16) 



Pdiff(x,y,z,x,y,z) 

= ((Ay-Ar (^Tr [ly^(x) W^^t(y) 



A^' 



:Tr 
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2 



X 



Ar2 



— Tr [W^\5c) W^if) 



W^,iz) - 1 



(A.17) 



Here we use the notations of (|4.4|J4.5| ) and F, M denote the analogous ex- 
pressions at transverse coordinates x, y. All terms quadratic in F (or F) are 
understood to contain the factor ln(xo/a;). Also we have not distinguished be- 
tween the true and the projectile biased inelastic amplitude, since it does not 
matter as long as we are interested in particles at rapidity far from the target. 



B. Application of Fierz identities 



In order to simplify the colour algebra in products of Wilson lines, we use the 
relation between adjoint and fundamental Wilson lines. 



2Tr 



fB.ll 



where the generators of the fundamental representation satisfy 



(B.2) 



We make extensive use of the Fierz identities 



Tij T^i - -^[^ii ^jk - J^^ij ^kl ] , 



(B.3) 
(B.4) 



The fundamental Wilson lines in (|B.4|) depend on the same transverse coor- 
dinate. This implies that transverse momentum integrations in cross sections 
of the type ( A.12| ) or ( A.13|) are needed for this Fierz identity to apply. The 
simplifications which we shall obtain in using ( [B.4|) are technically identical to 
those obtained in a Feynman diagrammatic language via Mueller's technique 
of cancellation between real and virtual (contact) scattering contributions. 
Indeed, in the Feynman diagrammatic language, the only prerequisite for ap- 
plication of Mueller's cancellations is again that parton lines are closed at the 
cut, i.e. that transverse coordinates of fundamental Wilson lines are equal in 
amplitude and complex conjugate amplitude^. The use of ( p.4|) will allow us 
to circumvent explicit diagrammatic calculations. 
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Using the Fierz identities ( |B.4| ), one can simplify e.g. the following contribu- 
tions to different DIS observables given in Appendix ^ 



Tr 



N' 



-Tr 



(B.5) 



Tr [Vr^(x) W^\y)] Tr [w^iy) W^\x) 
= 1 + W^,i^) W^.iy) , 



(B.6) 



Tr 



ly^^(x) W^{y)] Tr [t"^ Vr^(x) W^\y) 

W^,{z) W^,{y) 



W„i(z) 



2Ar 



2A^ 



+-W^,{z)W,j^{y)Wff{^) {d. 



if ace) {dbdf + ifbdf) 



(B.7) 



Tr 



W^a1(z) Tr 



4iV2 



+ 



SAT 
1 



iyi(z)W,'^(z)l^,^/y) (dace+tface) {ddbf + tfdb}) 



W^,{Z) W^IZ) W^Ay^) {daec + ifaec) (ddfb + ifdfb) 



+Yq ^^^^^^ W ^/^(y) 



X 



X 



^ae ^cg ~l~ {daem ~l~ if aem) idcgm ~l~ ifcgm) 



N 



Sfb Shd + {dfbn + iffbn) {dhdn + ifhdn) 



(B.8) 



C. Target averages 



To perform medium averages ((...)) over products of fundamental or adjoint 
Wilson lines, we expand the Wilson lines locally for each longitudinal position 
up to second order in the colour potential. For the last position ^ of the Wilson 
lines, one expands e.g. 

W^b{^) = Vl{z) (v-^^(e,z)/^^^-^ 5^5(^,0)) , (C.l) 
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W'\7,) = V' {z) ( 1 + ^ ^^(e, z) - ^ B{^, 0) ) . 



(C.2) 



Here, V"^ and denote the Wilson lines up to the last position ^. Also, we 
have anticipated that in the medium averages coulour fields are expanded only 
up to second order locally in the longitudinal direction. 



{{A'^{C,z)A\C,^')))^S^'BiC,z-z'). 



(C.3) 



In this way, one obtains e.g. for the average of two fundamental or adjoint 
Wilson lines 

((lTr[iy^t(^)iy%)])) 

(( llV [V^^i^) V^{y)])) (1 - m, 0) - BiC, X - y)]) 



exp[-CFv{x-y)] , 
((Tr [W^\y) W^{^)])), = exp [-C^v{^ - y)] 



(C.4) 
(C.5) 



Here, the last line is obtained by reexponentiation with the help of 

^(x) = I d^m,0)-B(^,^)] . (C.6) 



More complicated averages can be obtained in the same way: 

(( W^M^) Tr [W^\y) W^{y) T'] )) 
= ((KlWTV[y^t(y)r«T/^(y)r^])) 



X 



- ^ 0) + y - y) - y - x) - y - x)] 



-CF[Bi^,0)-Bi^,y-y)] 



_ -1 -%(„(x-y)+^(x-y))-(Cf-%-)^(y-y) 



(C.7) 



As soon as colour octets propagate along the products of Wilson lines, the 
calculations can become more involved. For the combinations of /-symbols 



^bd — £ £ 
•'b\d\ — Jbinbjdmd' 



(C.8) 



which appear in the averaged Wilson lines, it is useful to observe that both 
singlet and octet structures are eigenstates of propagation. 
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^bd f _ Ca r. 

■J bidi Jbdf — 2 Jbidif } 

nrbd J —^^ A 

bidi "'bdf - -TT "'bidif ■ 



(C.9) 
(C.IO) 

(C.ll) 



For an average of three Wilson lines, this leads to 

idace-^face) ( ( l^i(z) iy,^(z) W^,j(y) ) ) (d^f - tfbdf) 
= {dace - ifaee) (( Kt(z) ^ (z) V;^^(y))) 

3 



X 



Ca B{^, 0) 6b,b Sd,d 6f,f + B{^, z - z) 6f,f 



B{^, z - y) 6,,, + B{^, z - y) J^Zh K 



df 



{dbdf — ifbdf) 



{daee ' ifaee) (( KK^) K1 (^) ^^(y))) (4,../, - ifl 



bidifi) 

^5(e,0) + ^5(e,z-z) +^5(e,z-y) +^5(e,z-y) 



(iV^-l) (2iV--)exp 



Ca 

— Y b(z - z) + ^^(z - y) + ?^(z - y)] 



(C.12) 



More complicated is the average over the four adjoint Wilson lines given in 
(|B.8|) . Again, we expand the Wilson lines according to ( |C.1| ), 



((w^i(z)w^i(z)iye/(x)w^i(y))) 



(( K^(z) V^z) Fi(x) yj(y))) MfH 



(C.13) 



We introduce the shorthands 

Ra = -v{z - x) - v{z - y) , = -w(x - y) . 



(C.14) 



One eigenstate of M'^^j^ is rather easily checked, 

■^bdfh iffbndhdn ~l" dfhnfhdn) ~ Ra ifeandgen ~l" '^eara/gcn) • 



(C.15) 



To find the other eigenstates of M^^jj^, we introduce the vectors 

Ul = 6ea5ge , ^2 = ^ca^ge , % 
^1 feanfgen j ^2 feanfgen j ^3 

^ean^gen ; 



^ga^ee i 



fganfe. 



UJl dga^dgcn J U)2 df-njidg^n j ^3 dg^nd, 



ganj een 5 
gandecn ■ 



(C.16) 
(C.17) 
(C.18) 
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Colour algebra implies the identities 
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{U2 


- U-i) + W2 
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- Us) + Wi 


- W3 
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2 
N 




- Ui) + W2 


- Wi 



(C.19) 
(C.20) 
(C.21) 



In the six-dimensional space, given by the vectors (mi, ^2, M3, Wi, W2, W3), one 
finds by explicit calculation that the matrix M'^^^l acts hke (using Rd = Rb — 

Ra) 
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-Rd (I - f 

f {3Ra + Rb) J 
(C.22) 



To calculate the average in ([B.8|) , we have to propagate the vector 



■J^ Sfb Shd + dfbndhdn — ffbnfhdn 



— {Ui -U2+ Us) + Wi - W2 + Ws . 



(C.23) 



Finding the eigenvectors and eigenvalues of ( P.22|) with the help of Mathe- 
matica, decomposing the vector ( |C.23| ) in terms of the eigenvectors of ( p.22| ) 
and reexponentiating the obtained expressions, we find 



16 
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iV2 



e^^". (C.24) 
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This enters the calculation of the probability Pjf^^'^'^^'' in ( |4.15| ) and the diffrac- 
tive evolution equation ( |5.9|) . 



Finally we present averaged expressions for the single differential quark cross 
sections, obtained from ( |A.15| , |X?T7| ) by taking y = y. 



qq 



dp 



(C.25) 



dp 



l^~%l(^{x-y)+«;(x-y))-(cf-^) ^(y-y) ^ 



(C.26) 
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